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Abstract:  A  new  formulation  using  the  reduced  magnetic  vector  potential  and  the  electric  scalar  potential 
is  proposed  for  modeling  of  eddy  current  (EC)  problems  using  finite  element  analysis.  The  major 
advantage  of  the  proposed  formulation  is  that  it  eliminates  the  need  for  meshing  the  current  source.  In  EC 
nondestructive  testing  the  probe  coil  scans  across  the  sample  which  involves  motion  of  the  coil. 
Simulation  of  EC  scans  in  conventional  finite  element  models  requires  considerable  re-meshing.  In  the 
proposed  formulation,  the  mesh  and  system  matrix  remain  unchanged  with  different  coil  positions.  The 
incomplete  factorization  of  the  system  matrix  as  pre-conditioner  is  performed  only  once  during  the  entire 
scan,  which  results  in  a  significant  reduction  of  computation  time. 

Keywords:  Reduced  Magnetic  Vector  Potential,  Eddy  Current,  Finite  Element  Method 


1.  Introduction 

Finite  element  methods  based  on  potentials  have  been  successful  in  modeling  eddy  current  (EC) 
problems.  These  formulations  in  general,  require  the  discretization  of  the  current  source.  However,  in  EC 
nondestructive  testing  (NDT)  applications,  the  meshing  of  the  EC  probe  excitation  coil  becomes  very 
cumbersome  particularly  if  the  scanning  motion  of  the  probe  coil  is  to  be  taken  into  account.  A  simple 
example  is  a  coil  above  a  conducting  plate  with  a  circular  hole.  If  the  coil  and  the  hole  are  not  coaxial, 
meshing  becomes  complex.  As  the  coil  scans  the  sample,  re-meshing  the  coil  at  each  scan  position  is 
necessary  and  very  difficult.  It  has  been  shown  that  if  the  current  source  is  not  discretized  accurately,  the 
convergence  of  solution  becomes  slow  [1].  Hence  formulations  that  do  not  require  meshing  the  current 
source  are  highly  desirable.  Such  formulations  offer  many  advantages  as  listed  below: 

a)  Mesh  generation  with  complex  coil  shape  becomes  much  easier  because  the  current 
source  is  not  meshed  and  is  not  integrated  with  other  structures  in  the  solution  domain. 
The  coil  is  discretized  independently  for  numerical  integration  of  analytical 
expressions  of  the  field  due  to  the  source  current. 

b)  Convergence  of  the  iterative  solver  is  faster  because  the  current  source  can  be 
modeled  accurately. 
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c)  Solution  domain  can  be  reduced  because  the  secondary  field  generated  by  induced 
currents  and/or  magnetization  is  much  smaller  than  the  source  field. 

d)  Large  coils,  such  as  the  saddle  coil  used  in  the  magnetic  resonance  imaging,  do  not 
need  to  be  meshed  with  the  target  structure,  which  again  reduces  the  solution  domain. 

e)  It  is  easy  to  continuously  vary  coil  parameters,  such  as  liftoff  and  tilt  angle,  for 
parametric  studies  without  changing  mesh. 

f)  Scanning  of  coil  over  a  sample  involved  in  a  typical  EC  NDT  process  can  be  modeled 
without  re-meshing.  The  computation  error  due  to  re-meshing  is  removed. 

g)  Using  appropriate  formulation,  as  proposed  in  this  paper,  incomplete  factorization  of 
the  system  matrix,  which  is  the  most  time-consuming  process  in  the  solution  phase,  is 
performed  only  once. 

The  A,V-if/,(j)  formulation  employing  the  magnetic  vector  potential  A,  the  electric  scalar  potential  V, 
the  total  magnetic  scalar  potential  yr,  and  the  reduced  magnetic  scalar  potential  (j)  does  not  require 
meshing  the  current  source.  Unfortunately,  this  formulation  is  complex  and  more  importantly  it  gives  ill- 
conditioned  system  matrix  in  the  presence  of  high-aspect-ratio  elements  and  discontinuities  in  material 
properties  [2].  A  formulation  using  the  magnetic  vector  potential  A,  the  reduced  magnetic  vector  potential 
Ar ,  and  the  electric  scalar  potential  V  is  presented  in  [2] -[3].  In  this  paper,  we  propose  a  new  formulation 
using  only  the  reduced  magnetic  vector  potential  and  the  electric  scalar  potential.  With  this  AnV-Ar 
formulation,  the  system  matrix  remains  the  same  despite  changes  in  coil  position.  The  incomplete 
factorization  of  the  system  matrix  is  performed  only  once,  which  results  in  significant  savings  in 
computation  time. 


2.  Formulation 


2. 1  Magnetic  Vector  Potentials 

The  magnetic  flux  density  B  can  be  decomposed  into  two  parts,  namely, 

B  =  BS  +Br=jU0Hs+Br  (1) 

where  jU0  is  the  permeability  of  air;  Bs  and  Hs  are  the  flux  density  and  field  intensity,  respectively  due  to 
the  excitation  current  density  Js  in  free  space  and  Br  is  the  flux  density  due  to  induced  currents  and/or 
magnetization.  It  is  important  to  note  that  these  fields  satisfy  the  conditions  V  Bs  =  0  and  V  x  Hs  =  Js. 
Since  V  •  Br  =  V  B  -  V  •  Bs  =  0,  Br  can  be  written  as  curl  of  a  vector,  namely,  Br  =  Vx  Ar.  Ar  is  called 
the  reduced  magnetic  vector  potential.  We  can  then  express  B  in  terms  of  Ar  as 

B  -  HqH s  +  V  xAr  (2) 

where  Hs  is  calculated  using  Biot-Savart’s  law: 

LJ^xr~^dQ-  (3) 

4  7rJ^  \r-r\ 


In  (3),  Q'  is  the  volume  of  the  current  source;  r  and  r'  denote  the  coordinates  of  observation  and  source 
points,  respectively.  The  commonly  used  magnetic  vector  potential  A  is  now  called  the  total  magnetic 
vector  potential.  A  can  be  decomposed  as  A  =  Ar  +  As ,  where  As  is  the  magnetic  vector  potential  due  to 
the  source  current  in  free  space  and  is  expressed  as: 
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(4) 


The  A,V-Ar  formulation  is  presented  in  [2]-[3],  in  which  A  and  V  are  used  in  conductors;  Ar  is  used  in 
free  space.  The  two  vector  potentials  are  then  coupled  on  conductor  surfaces  using  As.  In  the  next 
subsection,  we  propose  the  new  AnV-Ar  formulation  that  does  not  use  the  total  magnetic  vector  potential. 

2.2  New  Ar, ,  V-Ar  Formulation 

The  new  formulation  stems  from  the  classical  A,  V-A  formulation  [4]: 

V  x  vV  x  A  +  jcooA  +  oVV  =  JS  (5) 

V  ■(jco<jA  +  oVV)  =  0  (6) 


where  v,  cr,  and  co  are  the  reluctivity,  conductivity,  and  angular  frequency,  respectively.  Equations  (5)-(6) 
are  used  in  the  whole  solution  domain;  whereas  only  (5)  neglecting  the  eddy  current  term  is  used  in 
nonconducting  region.  Substituting  the  decomposition  of  A  and  the  representation  of  Js  by  Hs  into  (5)-(6), 
we  have 


VxvVx^  +  jcooAr  -hoVV  =  Vx  Hs  -Vxi/Vx  As  -  jcooAs  (!) 

V  •  (-  jcooAr  -  oV  V)  =  V  •  jcooAs  (8) 

Applying  the  Coulomb  gauge  on^4r,  i.e.,  V  •  Ar  —  0,  and  assuming  constant  v  within  each  element,  (7) 
can  be  rewritten  as 


-  vV2Ar  +  jcooAr  +  oV  V  -  (l  -  vr  )V  x  Hs  -  jcooAs  (9) 

where  vr  is  the  reciprocal  of  the  relative  permeability.  Equations  (9)  and  (8)  are  the  governing  equations 
of  the  AnV-Ar  formulation. 

Note  that  the  weak  form  of  (8)  is  written  as 

f  (w,  •  jcooA r  +VV,  -(NV)dn  +  i  Nt  [-  jaxr{Ar  +  As )  -  oV  v]  ■  ndT  =  -  f  VV,  •  jcooA s dO. .  (10) 

J  Q  JF  J  Q 

By  neglecting  the  surface  integral  in  (10),  we  implicitly  set  the  normal  component  of  induced  currents  on 
conductor  surface  to  0. 

In  contrast  to  the  A,V-Ar  formulation,  the  advantage  of  the  new  formulation  is  that  the  total  magnetic 
vector  potential  is  not  used.  Therefore  there  is  no  need  to  account  for  the  coupling  of  A  and  Ar  on 
conductor  surface.  With  the  new  formulation,  the  system  matrix  is  the  same  for  different  coil  positions. 
Hence  the  incomplete  factorization  of  the  system  matrix  as  pre-conditioner  is  conducted  only  once  during 
the  entire  scanning  process.  The  cost  of  the  new  formulation  is  that  besides  Hs ,  V  x  Hs  in  presence  of 
magnetic  materials  and  As  must  also  be  evaluated.  However,  evaluating  these  analytical  values  takes 
much  less  time  than  that  required  for  performing  the  incomplete  factorization. 

3.  Results 

The  sample  and  defect  geometry  used  in  the  simulation  studies  is  illustrated  in  Fig.  1.  The  conducting 
plate  has  dimensions  of  72  mm  x  72  mm  x  4  mm  and  conductivity  of  2.741  x  107  S/m.  A  through  hole  of 
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radius  4  mm  is  placed  at  the  center.  A  crack  of  dimensions  4  mm  x  1 .4  mm  x  2  mm  is  located  either  at  the 
top  or  bottom  of  the  plate.  The  coil  having  100  turns  and  excited  at  3  kHz  has  inner  and  outer  radii  of  4 
mm  and  6  mm,  respectively  and  height  of  2  mm.  The  liftoff  distance  between  the  bottom  of  the  coil  and 
the  top  of  the  plate  is  1  mm.  The  coil  scans  across  the  hole  in  the  x  direction  (the  y  coordinate  of  the  coil 
center  is  always  0)  with  1  mm  step  size  parallel  to  the  plate  surface.  The  mesh  in  the  x-y  plane  is  shown  in 
Fig.  2. 


Scan  direction 

■  ■  ^ 


(a)  (b) 

Fig.  1.  Geometry  investigated,  (a)  top  view  (b)  side  view. 


-20  -10  0  10  20  -20  -10  0  10  20 
(a)  (b) 

Fig.  2.  Mesh  in  the  x-y  plane  (units  in  mm),  (a)  in  the  plate  with  hole  and  crack  (b)  in  air  with  coil. 


The  model  is  validated  by  comparing  the  model  prediction  of  impedance  of  coil  above  the  plate 
without  hole  and  crack  with  the  corresponding  analytical  value  [5].  The  analytical  solution  for  coil 
impedance  is  0.1641  +  j  1.83 19  Q  and  the  numerical  value  predicted  by  the  model  is  S0  =  0.1648  +  j  1.7 123 
Q.  This  demonstrates  the  validity  of  the  new  formulation. 

The  complex  EC  probe  signal  (coil  impedance)  as  the  probe  scans  the  sample  is  calculated  for  the  3 
cases:  S\  -  sample  with  hole  but  no  crack,  S2  -  sample  with  hole  and  surface  crack,  and  S3  -  sample  with 
hole  and  subsurface  crack.  The  changes  in  impedance  due  to  the  hole  and  crack,  i.e.,  Si  -  So,  S2  -  So,  and 
S3  -  S0  are  plotted  in  Fig.  3.  The  effect  of  the  surface  crack  is  evidenced  by  the  asymmetry  of  the  signal. 
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Subsurface  crack  is  also  detectable  in  the  phase  plot.  Computation  results  obtained  using  commercial 
software  VIC-3D®  are  also  shown  in  Fig.  3.  Our  results  match  well  with  VIC-3D®  results. 


(a)  (b) 

Fig.  3.  Changes  in  EC  probe  impedance  due  to  the  hole  and  crack,  (a)  magnitude  (b)  phase. 

New  formulation:  -  hole; - hole  and  surface  crack; - hole  and  subsurface  crack. 

VIC-3D®:  °  hole;  x  hole  and  surface  crack  ;  0  hole  and  subsurface  crack. 


4.  Conclusions 

A  new  formulation  using  the  reduced  magnetic  vector  potential  and  the  electric  scalar  potential  for 
modeling  EC  problems  is  presented.  Application  of  the  formulation  in  NDT  with  coil  scanning  motion  is 
studied.  The  major  advantage  of  the  formulation  is  in  mesh  generation  since  the  coil  is  not  required  to  be 
meshed  at  each  position  during  the  scan.  Furthermore,  the  system  matrix  remains  unchanged  during  the 
scanning  process.  This  results  in  a  significant  reduction  in  computation  time  since  the  time-consuming 
incomplete  factorization  of  the  system  matrix  is  performed  only  once.  A  more  extensive  testing  of  the 
model  on  multilayer  structure  is  currently  underway. 
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